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Abstract 

We propose an experiment involving a gaussian laser tunneling through a twin barrier dielectric 
structure. Of particular interest are the conditions upon the incident angle for resonance to occur. 
We provide some numerical calculations for a particular choice of laser wave length and dielectric 
refractive index which confirm our expectations. 

I. INTRODUCTION 

In previous articles [1,2], we pointed out the analogies between laser interactions with a dielectric block 
and non-relativistic quantum mechanics (NRQM) . Tests of the latter physics should be easier through 
use of this optical analog, if only because we do not have to deal with probability amplitudes. One 
example is the transition in NRQM barrier diffusion from complete coherence to complete incoherence. 
Resonance phenomena occurs for the former limit but is absent in the latter. For example, the 
condition for coherence in one- dimensional NRQM diffusion is that the spatial size of the single 
■ incoming wave packet be much larger than the barrier size [3,4]. With coherence a single wave packet 

^ f is transmitted. With complete incoherence, unlimited numbers of wave packets are transmitted at 

regular intervals. Whence, incoherence phenomena are absent when considering only a single incoming 
plane wave because of its infinite size. Within the laser/dielectric analogy the wave packet size is 
substituted by the laser beam transverse spreads. Interacting with a dielectric block, as expected 
from geometric optics, multiple reflections occur [5]. Coherence attains only if overlaps dominate 
in the outgoing laser spots. One important difference between optics and NRQM is that the optical 
situation is stationary. Time in NRQM is essentially replaced by a spatial variable in optics. This offers 
some significant advantages, e.g. even low intensity beams can be detected by long time exposures in 
detectors. 

The fact that optical experiments provide an easy way to test quantum mechanical predictions, 
since the early 1970's, stimulated experimental and theoretical studies on quantum-optical analogies 
[6-14]. For clear and detailed reviews on the correspondence between optics and quantum mechanics, 
we refer the reader to the works of Dragoman [15] and Longhi [16]. 

Since any realistic setup is three dimensional, we must also work theoretically in three dimensions. 
Even when all interesting physics is dominantly bi-dimensional, e.g. the y-z plane defined by the 
incoming central laser direction (chosen by us as the z direction) and the normal to a structured 
dielectric system. There is a benefit in the extension beyond a simple one dimensional study to 
more dimensions. It means that the variation in "energy" in one-dimensional NRQM [17] can now 
be achieved by merely changing the incident angle. Indeed, it is the normal momentum component 
which now plays the role of energy in one dimension. For example, in optics, we can, under certain 
conditions, transit from diffusion to tunneling by simply rotating the dielectric structure relative to 
the incoming beam. 

A dielectric-air-dielectric system mimics a potential barrier. Experiments realized by a double 
prism arrangement [18, 19] have recently discussed the phenomenon of frustrated total internal re- 
flection. In this paper, we go beyond a dielectric-air-dielectric system and consider a double prism 
plus a narrow slab separated by small air gaps of equal size (see Fig. 1). This reproduces a twin 
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barrier situation [20-22] which should exhibit one of the most intriguing rcsuhs of NRQM, the res- 
onant tunnehng effect [23,24]. Notwithstanding the expected attenuation of each barrier (air gap) 
upon the laser beam, for certain angular incidences the transition probability through twin barriers 
becomes unity, i.e. total transmission occurs. We determine the conditions for this phenomenon and 
reproduce sample results from a numerical calculation which necessarily takes into account the finite 
angular spread of any laser beam. By eliminating the intermediate slab, we reproduce the result of 
tunnelling through a single barrier [25,26]. The twin-barrier transmission results are attenuated by 
the entry/exit transmission amplitudes. All three must be folded into the laser momentum amplitude 
before integrating. This yields the outgoing laser beam amplitude and consequent intensity. For our 
analysis it is advisable to avoid the complication of multiple outgoing beams produced by incoherences 
within the dielectric prisms. This is probably best achieved by considering only small incident angles 
a. In principle it could also be achieved by reducing the depth of the prisms. This would facilitate 
the search for resonance at larger incident angles. Another effect that would also help in the use of 
larger incident angles, without lose of coherence within the prisms, is the fact that the y-profile of the 
laser beam is increased with a as 1/cosa. The x-profile is constant throughout and plays no role in 
the question of coherence. In the next section, we describe our theoretical treatment of the gaussian 
laser. In the subsequent section, we describe our proposed twin barriers structure. In section IV, we 
provide the transmission/reflection probabilities for entry, exit and for the twin barriers. Section V 
contains our numerical transmission results for a chosen laser wave length and dielectric refractive 
index. The results confirm our phenomenological expectations. Our conclusions are drawn in the final 
section. 



II. THE GAUSSIAN LASER 

Discussion of the laser is simplified by first concentrating upon the electric amplitude it produces and, 
wherever possible treating the amplitude as if it where a scalar. This amplitude is a convolution in 
wave number of plane waves. For different plane wave directions the polarization vector necessarily 
changes since it lies in the plane perpendicular to the wave direction k. However, for a sharp gaussian 
in kx and ky we may assume to a good approximation that this vector (if linearly polarized) points 
in a fixed direction, 

E{x,y,z) = {Es{x,y,z) , Ep{x,y,z)} , 

where Es is the component of the electric field perpendicular to the plane of incidence (s-polarization) 
and Ep is the component perpendicular to this plane (p-polarization) . 

In this paper, we shall present the discussion of resonant tunneling for a gaussian He-Ne laser, 
A = 633 nm, of width Wo(= 2mm) which, in free space, propagates with principal axis along the 
^-direction, i.e. 



G{kx,ky) = exp 



^ ( 
4 V 



(1) 



To simplify our presentation, we shall suppose that the laser beam passes through a linear polarizer 
which selects the s-polarization. This allows, as shown in the next section, to directly use the QM 
results. The electric field can be thus treated as a scalar, 



E{x, y,z) = { Es{x, y, z) , Ep{x, y,z)} = { E{x, y,z),0} 



where 



E{x,y,z) 



47r 



dkx dky G{kx ,ky)e' . 



(2) 



Observing that fcwo = 27rwo/A 3> 1, we find the following analytic analytic expression for the free 
space intensity. 



I(a;,2/,z) = lo 



w^(z) 



exp 



X +y 
w^{z) 



Kz) 




(3) 
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III. LAYOUT OF TWIN BARRIER APPARATUS 



In Fig. 2, we show the y-z overview of a prism which forms a part of a twin barrier structure. The 
incident ray corresponding to the central direction of the laser travels along the z-axis. The incident 
angle with the normal, z*, of the leading face of the dielectric is denoted by a. The internal direction 
P satisfies Snell's law, 

sin a = n sin j3 . (4) 

Before to derive the conditions for total internal refiection (TIR), let us establish the connection 
between optical and quantum mechanical problems. This connection will be the starting point for the 
calculation of the transmission probabilities presented in the next section. 

The Fresnel formulas for the reflected and transmitted waves obtained when s-polarized light moves 
from air to a dielectric of a given refractive index n are [5] 

cos a — ri cos (3 2 cos f3 . 

Rs = ■ ^ and Ts = ■ . (5) 

cos a + n cos p cos a + n cos p 

The air dielectric system is, for s-polarized incoming waves, the NRQM analog of a step potential. 
Indeed, for a step potential with a discontinuity perpendicular to the z^-axis we obtain the following 
reflection and transmission coefficients [17] 

Rstep = ^ and Tstep = T ~ — ■ (6) 

Now, observing that for plane wave 



^z, = k cos a and = yri'k^ — ky^ ~ \J nk^ — k^ sin^ a = nk cos /3 , 

we find Rstep = Rs and Tstep ~ Ts- This analogy can be claimed even if there are subsequent 
structures for we now know that all phenomena can be calculated as successive events even when 
coherence reigns. For the above situation there are always reflected and transmitted beams since the 
momentum vector is greater within the dielectric than in air. We are in the analogy of QM diffusion. 

Passing from the dielectric to air, we have a situation where for angles greater than a critical 
incident angle 7c, 

n sin 7c = 1 , (7) 

total internal reflection (TIR) occurs because Snell's law cannot be satisfied. If the air region is of 
finite dimensions, we are in the analogy of QM tunneling. 

The twin barrier structure, see Fig. 1, is created by air zones separating two prisms from a nar- 
row slab made of the same material as the prisms. Since the resonance phenomena that we seek to 
reproduce requires coherence, the total distance between the two prisms must be much smaller than 
the laser beam size (with one notable exception - head on incidence). For certain angles, as we shall 
demonstrate in the next section, transmission through the twin barriers will be unity (resonance). 
However, if the barriers are not of the same size resonance is attenuated. The direction normal to 
this twin barrier structure we call z. The overall transmission (refiection) amplitudes comprising of 
the entry step, twin barrier and exit step must multiply the gaussian modulation before integrating 
over wave number plain waves. In general this can only be done numerically. Obviously, the entry 
and exit step amplitudes will produce a reduction in laser intensity even at resonance. This should 
not be a serious problem for a steady state intensity, which can be measured over an arbitrary time 
interval. As a check of our results, the simpler single barrier amplitude can be derived from the twin 
barrier one by setting the intermediate-slab thickness to zero. The resulting amplitude is that for a 
barrier width equal to the sum of our twin barrier sizes. 

Before we calculate the transmission amplitudes, it is useful to derive the conditions for TIR. To 
do so we must consider two distinct cases. 
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• Case (a). Shown in Fig 2a. From the internal triangle with angles 7, 7r/2 — /3, and 7r/4, we find 

7 = 7r/4 + /? . (8) 

TIR occurs for values 7 > 7c. This implies the following constraint on /? 

sin2/3 > sin(2 7c - 7r/2) . (9) 

Observing that sin 2/3 > 0, the previous condition is always satisfied when 7c < 7r/4 since the right 
hand side is negative. Thus for n > \/2 only TIR occurs. Let us now investigate the constraint on a 
for TIR in the case n < \f2. From Eq.(9), we get 



sin2/? > 2 sin'7c - 1 = cos2/3<— ^ sin^ a > . (10) 

■n? 2 

Summarizing, in case (a), the conditions on the incident angle a which guarantee TIR are 

n>V2: Va , 

^ . ln'-2V^^ir^ (11) 
n < v2 : a > arcsm y . 

• Case (b). Shown in Fig 2b. This is obtainable by flipping the prism (and consequently the full 
dielectric set up) with respect to the incoming central laser direction. Alternatively, and perhaps 
more practically, it is obtained from the previous case by flipping the laser beam about the normal to 
the prism. From the internal triangle with angles tt/2 — 9, tt/2 — /3, and 7r/4, we find 

7:=7r/4-/?. (12) 

TIR (7 > 7c) imposes the following constraint on /3 

sin2/3 < sin(7r/2 - 2 7c) . (13) 

Observing that sin 2/3 > 0, the previous condition is never satisfied if 7c > 7r/4 since the right hand 
side is negative. Thus for n < \/2 only diffusion occurs. For n > \/2, from Eq.(13), we have 



sin2/3<l-2sin^(?c = ^^^ => cos2/3>^^^^^ sin^ a < ^ "^"^ V (14) 

Finally, in case (b), the conditions on the incident angle a which guarantee TIR are 

n<V2 : , 

/- . /,„ -2V^^^^ (15) 

n > V2 : a < arcsm ■ 



IV. TRANSMISSION PROBABILITIES 



The principal transmission amplitude emerging from the far edge of the dielectric structure as sketched 
in Fig. 1, is a product of three separate transition amplitudes, the entry and exit transitions and the 
twin barrier transition (assumed coherent in this study). Reflected amplitudes also occur at each 
stage and will in general give rise to secondary and higher transition amplitudes which we ignore for 
simplicity. Since the perpendicular momentum or wave number is of such relevance much of what 
follows involves the use of appropriate rotation matrices. We start with the entry step and consider, 
prior to convolution, an incoming plane wave. 

• Entry [air /dielectric interface]. 

At this interface the standard one-dimensional step results hold with the relevant momentum given 
by the z* component of k. The direction is normal to the interface and in terms of our original 
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axis (chosen so that the principal direction coincides with the z-axis) it corresponds to a rotation of 
a around the x-axis, 

' z* \ _ / cos a — sin a \ f z 
y* ) \ sin a cos a J \ U 

The wave number components of the incoming laser beam which propagates in air can be expressed 
in the new axis in terms of the wave number fc, 

^y. = ^2 sin a + ky cos a — k sin a + ky cos a + 0[k^, ky] , (16) 
= kz cos a — ky sin a = k cos a ^ ky sin a + 0[k^, ky] . 

The wave number components of the laser beam which propagates in the dielectric are 

Qx^ — kx^ ; 

<ly, = ky*j (17) 

Qz, = yj/n^k^ — kl^ — k^^ = nk cos (3 — ky tan/3 cosa + 0[kl,ky] . 

Now the expressions for the transmission (reflection) amplitudes have phases depend upon the value 
of z* for the step. Our laser formulas are written with the origin as the central point of minimum 
transverse size of the laser beam. We denote the z* distance of the interface from this by a* (see 
Fig. 1). With these definitions, we recall the step results 

kz, + qz, * k,, + qz, 

• Exit [dielectric/air interface]. 



For the exiting stage dielectric/air interface at z* = 6*, we need only interchange in the entry results 
kz, with qz, and change a* with 6*, 

= e^'-?-^- and T^' = ^H,.,-k.,)b, (^g) 

Qz, + kz, qz, + /c.. ^ ' 

For our calculation of the transmission amplitude, we shall need the product of the entry and exit 
transmission amplitudes. Tee, 

TEE(fc.,fc.) = <-"'t::^' = ^i,.,-^.,)D, (20) 

(fcz. + qz,) 

where = 6* — a*. In Fig. 3, we show the values of the above product for different values of refractive 
index n as a function of the incident angle a. The values of k^ and ky have been chosen as zero here, 
so these curves strictly apply to the central laser ray. The dependence is mild (particularly for small 
incident angles) compared, for example, to that induced by resonance (see below). Its main effect 
will be that of reducing the emerging laser beam intensity compared to the ingoing laser beam, even 
at the peak of resonance when the twin barrier amplitude is unity. As pointed out in the previous 
section, this should not constitute a serious problem. 

• Central structure [dielectric/air /dielectric/air /dielectric interface]. 

For resonance phenomena in tunnelling, we consider the central structure as small or comparable to 
the laser dimensions and give the overall transmission (reflection) amplitudes for this. Before doing 
so it is convenient to rotate the axis anew about x to define the z direction normal to the central 
structure 

( " 71 ( 1 1 ) ( y* ) ■ 
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In terms of the new axis, the dielectric wave numbers become 



Qy ~ (gz^ + Qy, )/-\/2 = nfc sin7 + fcy COS7 cosa/cos/3 + 0[fcj, fc^] , 
qs = {qz, — qy,)/ \/2 ^ n k cos J ~ ky siw-f cos a / cos 13 + 0[kl, ky] , 



(21) 



and the air wave numbers 



qx ■ 

q ~y . 



) ni '2 ~ 1 R 2.- 2 . sm7COS7 cosg 

ki = Jk^ ^ q^- q- = ky'l- sm 7-71 fc.y —^====== ^^^-^ + 0[k^, ky 



(22) 



y/l - sin' 7 COS/? 

For n sin 7 > 1, we have tunnchng. In this case, we can rewrite kz = i\kz\ where 



1; I ; / , • 2 7 , , sin 7 COS 7 cos a 

\kz I = k Vn- sm 7 - 1 + n ky : — — 

sm 7 — 1 cos p 

The single barrier reflection and transmission coefficients arc [17] 
R- 

where 



^y] 



j^l±J^ sinhfl/c^lL) e"^/ J- and T = V -7" 



(23) 



ql + M 

2q-z\ki\ 



sinh ( L 



(p = arctan 



2qi\k-z\ 



tanh [\kz\L 



and L is the dimension along the z axis of the air gap (see the figure below) . The amplitude due to 
a back and forth reflection within the dielectric slab region produces the following loop factor 

2 ~ 2 ~ 

i? exp[2 i 92 5] X i? exp[ — 2 i a] = i? exp[ 2 i gj(6 — a) ] = — exp[ 2 i (gj D + (^) ] , 

where D is the dimension along the z axis of the dielectric slab (see amplification in Fig. 1). Thus 
performing the sum over successive loop contributions to the twin barrier amplitude, Ttb, yields 



Ttb [kx , ky 



T 



l + \R\' exp[2i(gjL» + ip)] 



(24) 



The above result can also be calculated through a standard matrix approach. A check of this result 
can be made by setting D = {). This yields the single barrier result for a barrier size of 2 L to be 
compared with that for the size L barrier results used in the above derivation. The results are in 
accord. 

The resonance conditions are achieved when the phase factor multiplying \R\ equals —1. The 
denominator is then equal to |r| and the twin barrier amplitude reduces to a mere phase and hence of 
unit modulus. In general the denominator wilHar exceed the numerator, which for sufficiently large L 
is proportional to a decreasing exponential in L. In Fig. 4, we show some explicit examples of this phe- 
nomena. In this figure, we have set arbitrarily fc^, = and ky = and plotted ITtbI versus the incident 
angle a for a refractive index n = Note that this angle is not the impinging angle upon the twin 
barrier structure itself (see Fig. 2). The three curves are for {kL, kD} = {1.5, 100}, {1, 100}, {1, 500}. 
In each the resonance structure is clearly manifest. The peaks are very close to each other in the 
latter case. The first two plots demonstrate that for a fixed D the resonance peaks are narrower as L 
increases. As an aside, we recall that resonance requires the existence of a non null D. The resonance 
condition cannot be obtained by alone. For tunnel resonance one needs multiple barriers. 



6 



V. NUMERICAL ANALYSIS 



Our interest in this paper is the outgoing transmission amphtude from the twin barrier dielectric 
set-up described in the previous section and sketched in Fig.l. In particular, we are interested in 
the resonance phenomena displayed within the tunnelling regime, i.e. for TIR. TIR introduced for a 
step potential is somewhat of a misnomer for finite barriers. Indeed resonance, when it occurs, is the 
extreme counterexample since it represents total transmission not total reflection. The electric field 
for the transmitted laser is obtained by multiplying the gaussian G{k,j.,ky) in Eq.(2) by T^^ik^Tky) 
and TrY-aiJ^xi f^y)j 



ET{x,y,z) 



w= [ 
47r J 



pi (k^ x+ky y+k^ z) 



Observing that fcwo ^ 1 and that the dependence on kx in Tee ?tb is only of second order, whereas 
the ky dependence is of the first order, we can simplify the transmitted electric field as follows 



\Er{x,y,z)/E,\ 



2 y/T:w{z) 



exp[— x'' /w^{z)] 



dky Griky) exp 



kyV-^^ 



(25) 



where 

Griky)=eM~kywl/A] rEE(0,fcj,)T™(0,fc.y) . 

In Fig. 5, we plot |Gr(fcj,)| as a function of ky around the angular resonance for the following choices 
{kL.kD} = {3,10^}, {1,10^}, {1,10*}, {1,5 • lO"}. Also plotted in the same figure is the incoming 
gaussian with k^ = 0, i.e. G(0, ky). The effect of rEE(0, ky) TtbIO, ky) reduces the peak values of each 
curve and for {/cL, kD\ = {3, 10 } shifts the peak significantly to a positive value of ky. In the case 
with {kL,kD} = {1,5 • 10 }, we see the presence of multiple peaks. Observe that the values of L are 
very small if compared to the laser beam size, kL = 1 corresponds to L = 0.1 /zm. 

The spatial size of the laser is set by wq and is of 2 mm. For distances between laser source, 
dielectric structure and measurement apparatus of the order of a few meters we can readily neglect 
spatial spreading due to the gaussian, although this does not necessarily apply (as we shall see) to 
the effect of twin barrier resonance upon the beam. As we have argued in depth elsewhere, resonance 
phenomena requires a twin barrier structure smaller in size than the beam size. This means that 
Z? + 2L <C 2mm. Since for other reasons L must be much smaller (of the order of A) this limits D 
to less than a mm. For greater values of D the resonance peaks proliferate so that the gaussian laser 
would encompass multiple resonance peaks. In this case varying the incident angle produces little or 
no consequences and the resonance effect is averaged out. 

With the same selection of values of {L, D} used in Fig. 5, we have performed a numerical calcu- 
lation of \Ej,{x, y, z)\^. In Fig. 6, we plot the exit contour of the laser intensity in the x-y plane at the 
1/e^ width. The results are, of course, correlated to those of Fig. 5. For example, a narrow curve in 
ky such as in Fig. 5a produces a large spread in the y contour in Fig. 6a. The a;-spread is the same in 
all of the Fig. 6 because the x-integral in Eq. (25) factorizes (after ignoring second order corrections in 
kx)- In Fig. 6d, we display for comparison the incident laser beam with peak intensity Iq. The valures 
of the peak intensities (all close to x = y = 0) are also given in each figure. Fig. 5b corresponding to 
{kL, kD} = {1, 10 } is notable because it almost reproduces the incident laser reduced in intensity by 
a factor of 1/2. We do not display the more complex case corresponding to Fig. 5d with its multiple 
peaks. 



VI. CONCLUSIONS 



We have described in this paper is some detail an optical analogy of the twin barrier resonance 
experiment. The presentation has been made in three dimensions for an incoming gaussian laser 
beam characterized by a width of Wq = 2 mm. The laser wavelength has been set at A = 633 nm and 
the dielectric refractive index has been arbitrarily set at n = a/S. The twin barriers are reproduced by 
identical air intervals encompassing a narrow mm) dielectric slab. The width sizes which allow for 
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a clear angular separation of resonance peaks are of the order of several fj,m. Concentrating upon the 
first resonance peak, we find a very narrow angular spread comparable to that of the laser beam itself. 
To clearly see the resonance peak, our primary objective in this proposal, we must have a laser spread 
of the order or even smaller than the resonance peak. In the contrary case, the laser beam would 
encompass multiple peaks and varying the incident angle a would not result in significant differences 
in the outgoing beam. 

Once the resonance structure has been confirmed and assuming the numerical estimates of the 
previous sections are also convalidated, we can use the angular separation between resonance peaks to 
measure any of the parameters assumed in this work. For example the distance D or/and the barrier 
widths L. The values of the resonance angles may also be used to measure the laser wavelength 
if all the other parameters are known. In principle, it should also be possible to detect variation 
in the twin nature of the barriers. This could be caused both by modifying the equal sizes of the 
barriers or by modifying one or both of the refractive indices of the air barriers, e.g. by density, 
temperature or chemical nature of the "air" enclosed. Also probable is the loss of resonance due to 
external disturbances such as vibrations. While these are annoying characteristics for any experiment 
they also open the possibility of practical applications. 
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Figure 1: Schematic presentation of a double prism plus a narrow slab separated by air gaps. 
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Figure 2: Schematic presentation of a diagonally sliced dielectric prism. In the case (a), the condition 
n > \pl guarantees TIR for all the incidence angles (0 < sin a < 1). In the case (b), we have TIR for 

n > V2 and < sin a < \J n - 2\/n - 1/^2. 
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Figure 3: The modulus of the entry exit transmission amphtude for the central (s-polarized) laser ray, 
kx ~ ky = 0, plotted as a function of the incident angle, a, for different values of refractive index, n. 
The continuous line corresponds to the refractive index, n = -^/S, used in our numerical studies. For 
small incidence angles, the variation is very smooth and less than 10%. 
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Figure 4: The modulus of the twin barrier transmission amphtude for the central (s-poralized) laser 



ray, 



0, plotted as a function of the incident angle, a, for different values of the twin air 



gaps, L, and the dielectric slab, D. The resonance peaks are narrower as the air gaps increase and 
the distance between any two peaks decreases by increasing the dielectric slab dimension. 
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Figure 5: The incoming gaussian distribution, G{0,ky), and the transmitted distribution, \GT{ky)\, 
for s-porahzed laser beams of size Wq = 2 mm are plotted as a function of ky for different values of 
the twin air gaps, L, and the dielectric slab, D. The incident angle, a, is chosen around the first 
resonance for the first three cases, see plots (a,b,c). In the last plot (d), it is manifest the presence of 
multiple peaks. 
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Figure 6: The exit contour of the (s-porahzed) laser intensity at 1/e^ width is plotted in the x-y plane 
for different values of the twin air gaps, L, and the dielectric slab, D. The distance from the source 
of the laser, z-axis, and the dielectric block dimension, z,-axis, are respectively chosen to be 500 Wq 
(1 m) and 50 Wq (10 cm). The values of the peak intensities are given in terms of the incoming laser 
beam peak intensity Iq. 
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